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Introduction 

The objective of this report is to evaluate quantitatively the relative 
merit of three-component magnetic data over the total field measurements. At 
the present time, most magnetic surveys conducted for mapping subsurface 
geological structures are limited to measurements of the total magnetic field. 
This field is composed of the main field of the earth and the anomalous field 
cuased by crustal magnetized bodies. Since the anomalous field is normally 
much less in magnitude than the earth's main field, for all practical purposes 
It can be assumed to be measured in- the direction of the main field. This 
kind of observation vith total field magnetometer does not contain any infor- 
mation regarding change in direction of the anomalous field in space. This 
limitation of the data introduces ambiguities in the results of interpretation 
and does not permit a unique determination of any parameter of the causative 
bodies without additional geological and geophysical information. Recent 
theoretical studies (Bhattacharyya, 1967 ) Indicate ‘that if the three-component 
magnetic data are available, the direction and magnitude of the magnetization 
vector can be computed without any assumption on the geometry and orientation 
of the causative body. 

Measurements of one single component of the magnetic field anomaly over 
a horizontal plane contain Information about the other components. The 
possibility of deriving three-component data from one-component measurements 
was first suggested by Vestine and Davids (l945). Hughes and Pondrom {l9hl) 
showed that the vertical and horizontal anomalies can be obtained as surface 
integrals from the measurement of the total field anomaly. In practice, the 
development of a numerical method for obtaining three-component magnetic data 
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from measurements of one single magnetic component is based on transformation 
of the total -field data in the frequency domain. A method for expressing the 
data in terms of double Fourier series has been developed (Bhattacharyya, 19^5 )» 
With the help of this method, Lourenco (l972) and Lourenco and Morrison (1973) 
have computed the three-components of the magnetic field from the total field 
data obtained over a finite area. However, the accuracy of the calculated 
field components depends not only on the numerical method employed but also on 
the assumed direction of the geomagnetic field. It is very interesting to study 
the effectiveness of the calculated field components in the computation of 
the parameters of bodies causing magnetic anomalies. 

In the quantitative treatment of magnetic data, the first and second 
derivatives are usually evaluated. Second derivative maps are found to be 
very useful' in approximately delineating the boundary of the causative body. 

The objective of the present study is to evaluate the relative merit of three- 
component magnetic data over total-field meas'urements. Estimates are to be 
made of the. errors introduced in computing field components from total field 
data. The real test of the relative merit of the two sets of data is provided 
when they are used to determine the parameters of disturbing bodies. 

In chapter 1 the numerical method and results of obtaining three 
field components from total field measurements, using double Fourier series 
expansion, are presented. The horizontal first and second derivatives of each 
component are obtained by the method of bicubic spline interpolation and they 
are also discussed in the first chapter. The total field data for this study 
are computed over a finite area with a 6U x grid points, using a theoretical 
model of infinite prism. The field components and their derivatives, as 
calculated with the total field data, are compared with the exact values and 
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the shortcomings of the i^se of the total-field data are reported. 

In the second chapter of this report the expressions for moments of 
the anomalous field components over a finite area are given. The magnitude 
and direction of the magnetization vector are evaluated -with the help of 
moments of computed and exact three-component magnetic data for various 
prismatic bodies. The results of calculating the-magnitude and direction of 
the magnetization vector, using both computed and exact field component data, 
are uiscussed in this reuort. 

r- 

The study presented in this report shovs that the vertical component 
of the magnetic field calculated from total field observations is more accurate 
at higher geomagnetic latitudes than at lower latitudes. The opposite is the 
case for the horizontal components. The errors in determining the magnetization 
vector directions with the help of calculated field components are significantly 
large over most of the range of variation of declination and inclination of 
the vector. This observation demonstrates conclusively the practical limitations 
of computing field' components from total field data even under the best of 
conditions. The errors, as noted before, will definitely be larger for actual 


field observations. 
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CHAPTER 1 

MAGNETIC FIELD COMPONENTS AND THEIR DERIVATIVES 


1, Method for computing field components 

A method for computing three field components from the total magnetic 
field observations as described by Lourenco (l972) and Lourenco and Morrison 
( 1973 ) is presented in this section. Let us adopt a Cartesian coordinate 
system such that in the plane of measurement the x-axis points northward, the 
y-axis eastward and the z-axis vertically downward, as shown in Figure 1 . 

The magnitude of the anomalous magnetic field vector is generally very 
small compared to the geomagnetic field vector. This is why the direction of 
the resultant of the two vectors is practically invariant over the area of 
measiirement and is assumed to be the same -as that of the geomagnetic field. 

If AV represents the anomalous magnetic potential, at'.all points of free 
space it satisfies Laplace's equation: 

(1-1) '0'^ - 0 . 

The total intensity anomaly in the direction of the geomagnetic field is given 
by the negative derivative of the anomalous magnetic potential: 


(1-2) 
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where t is the unit vector in the direction of the geomagnetic field. 
The field components are given by the following equations: 
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^ere H , H and H are the field components along the x, y, and z directions 
x’ y z 

respectively. 

Since -both AY and satisfy Laplace .equation, the expression for 
the total field anomaly may be written in the following way in terms of a finite 


harmonic series expansion: 
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where L . L are the fundamental wavelengths along x and y axes respectively and 
X* y 

W„ is the discrete Fourier Transform of AT. 

Tmn 

In practice, the total magnetic field AT(x,y,o) is represented in terms 
of a grid of discrete values 
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in terms of a finite harmonic series expansion 


can also be represented analytically 
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respectively. 

Introducing (l-6) and (l-8) into the equations (l-3) , and Cl-5), 

we get the following relationships 
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where and F^, angular frequencies along x and y axes respectively, 
given by 
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where c<'» and ,-y' are the direction cosines of the geomagnetic field. 

It is to be noted that the expressions for W , W , and V nr(^ 

xmn’ ymn’ "*zmn 

undefined when m and n are simultaneously equal to zero. Therefore V ¥ 

xoo ’ yoo * 

^zoo calculated by other means as shown in the following 

paragraphs . 

The Fourier transform W(u,v) of a function of two variables F (x,y) is 
given by: ^ 

iv(u^\/)r JJ l^~ ITT i ( It jC ~h vy )~j d-7C ^ly ^ 

where u and v are the wavenumbers in the x and y directions respectively. 

From the above equation, we can obtain the following equation 
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It can be shown ^Bhattacharyya, I967) that 
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Similarly we have W^(0,0) = w (0,0) = w (o,0) = o. 

In practice, it is impossible to perform the integral over an iufinit 
plane and so the integral of the above equation will never be exactly equal 

to zero. Rewriting the equation in summation forms in terms of a finite grid 
Spacing, we have 
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Choosing the origin of the Cartesian coordinates at the lower left of 

the map and assuming that H k w a. 

® xjk’ yjk’ zjk ^ along the boundary, 

from equation (l-8) we get the following equation: ' 
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Thus and can be obtained from equations (l-8), (l- 9} and (l-lO) 

if the Fourier series expansion of the total field observations is available. 
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Eyaluati.on oS fi.&ld component-s 

In this section, the accuracy in calculating the field components from 
the total field data by the method described in the previous section will be 
examined with the help of a theoretical model. 

The model chosen is an infinite prismatic body as shown in Figure 1. 

The, dimensions and the declination and inclination of the magnetization vector 
of the prism are: 


2a = 16 , 2b = 8, h = 3, I = 20°, D = 0°, I = 60 ° D = lt5° 

o ’ O ’ 


where 2a and 2b are the horizontal dimensions along x and y axes respectively, 
h is the depth to the top of the body. I and D are respectively the inclination 
and declination of the geomagnetic field respectively. It is assumed that the 
prism is uniformly magnetised and the intensity of magnetization Is unity. 

The exact values of the field components of the prism are computed from 
the following e.iuationB (Bhattacharyya, 1964). 
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is the intensity of magnetization of the prism and c? 
direction cosines of the magnetization vector. 

The total field component is given by 
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where ^ y' are the direction-cosines of the geomagnetic field. 

In practice, the total field anomaly is represented as discrete samples 
over a finite area. The discrete Foiirier transform of these sample values is 
calculated ’vrith the help of the fast Fourier transform (FFT) computer algorithm 
(Cooley and Tukey, 1965). 

The anomalous field and its components were computed on a 6>h "by array 
of points with a unit grid interval. All of these computations were carried out 
with the help of a CDC 6ii00 computer in the Computer Center, University of 
California, Berkeley. 

Figures 2 and 3 show the exact and computed values respectively of the 

x-component of the field. The rectangle in broken lines represents the horizontal 
cross-section of the' prism. The large square contains 6h by Sh sample points. 

On the left hand side of each diagram are indicated the declination and inclina- 
tion of the geomagnetic field. Triangle represents the point of maximum value 
of the field. 

The anomalies in the two figures are practically Identical. The differences 
in the anomalies are noticeable when the field component is about one-half of its 
maximum value. The exact field attenuates slightly faster than the computed field. 

. Figures k and 5 show the exact and computed values respectively of the 
y-component of the field. Both the figures are practically identical. However, 
the computed field in this case falls off a little faster than the exact field 
and the discrepancies in the two figures become noticeable at small values of 
the field component. 

The characterisitics of Figures 6 and 7 showing the exact and computed 
values of the z-component are also the same as those in the above two figures. 

The disparity between the two figures can be seen only when the magnitude falls 
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off to about one-sixth of the juaxitnum value. 

figures 8 to 10 are the profiles along the K-S and E-W directions 
through the maximim response points. In each figure the exact and computed 
fields are marked hy syrobols (l) and (2) respectively. The horizontal 
axis of each figure represents the relative position on x and y directions. 

The vertical axis shows the magnitude of the field component along the profile. 
These profiles are useful for comparing the exact and computed fields in detail. 
As noted in the figures, the anomalies match remarkably well near the central 
parts where the magnitudes are high. However, in other regions there is 
noticeable discrepancy. For the particular model chosen, the match between 
the exact and computed values is the worst in the case of 

The prismatic ihodel in Figure 1 Is assumed to be located at a high geo- 
magnetic latitude, = 60° . If the model is placed at lower latitudes, it can 
be safely concluded that the differences between exact and computed z -component 
field will increase. However, because of the change in the characteristics of 
the anomaly with geomagnetic latitude, the accuracy in the computation of the 
horizontal components cannot be predicted with confidence. 

In order to study the effect of the variation of the geomagnetic field 
alone, all the parameters of the 'model in Figure 1 are kept fixed and only 
Iq and are changed to 30° and 15° respectively. Figures 11 - 13 present 

the profiles of exact and computed X, Y, Z field components along the N - S 
and E-W directions through the maximum response points. A comparison of 
Figures 8 - 10 and 11 - 13 indicates an increase in discrepancy between the 
exact and computed values of the field components. 

Figures l4 - l6 show the profiles of the field components for a model 
with the following parameters: 
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2a.= 8, 2b = 4 , tL = 3, I = 40*^, E = I = 30°, D =15° 

o o 

•V 

llie objective for the choice of this model has been to determine changes in 
the accuracy of the computed field components due to variation in the geomagnetic 
field as veil as in the horizontal dimensions of the prism. The effect of higher 
inclination of the magnetisation vector compensates to a certain degree that of 
the lover geomagnetic latitude in the present example. By comparing Figure 8 
vith Figure l4 and Figure 9 vith Figure 15 ve find that for x" y- components 
of the field, better results are obtained in the present case than in the first 
case., Hovever, for the z- component, the differences betveen exact and computed 
field values near the central part of the anomaly are higher in Figure l6 than 
in Figiore 10. 
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3. Horizontal Gradients of Field Components : 

In the previous secton, we discussed the method for obtaining the field 
components from total field measurements using finite harmonic series expansion. 
This section deals with the comparison of the horizontal gradients of the 
computed field components with the exact field gradients. Spatial filters were 
not used for this computation becasue the first horizontal derivative has an 
imaginary response in the frequency domain. 

Bhattacharyya (1969) used bicubic splines representing a given set of 
magnetic observations for calculating the horizontal derivatives. This method 
is used to compute the derivatives of components. A brief treatment of this 
method is given in the following paragraphs. 

If a function g (x,y) with two variables x and y is given at every mesh 
point of a grid system, such as 



then the piecewise bicubic spline function has the following form: 
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can be easily computed if the coefficients a ^ are available. 
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The parameters p, q. and s can "be evaluated hy the following steps: 

Ci) calculation of p. . 

ij 

For each of the lines -j ~ C: , i ^ x ^ - first derivatives 

p^^ should satisfy the equation. 


Cl-IT) 


h- . *f . 4- , 




3 

, M-l 3 

At boundaries j when i = 0 and i = M, the following condition holds: 


< = h ^ , ' 


Ci-l8) 


) = ° , 


^ Z Cy M ) 


f - I, X, ■■■ /^ • 


(ii) calculation of q. . 

The equations for calculating q are similar to those of p. The equation 
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(iii) Calculation of g • ^ 

For computing s. .» the values of p and q in equation (1-1T) aJ^id (l-19) 
are used. 
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The above steps are necessary'- for computing p, .q.j and s. The detailed 
treatment is presented by Bhattacharyya (1969)- 

With the help of the values of p, q, and s» it is easy to compute the 
coefficients of the piecewise bicubic splines. It is also very straightforvrard 
to calculate the derivatives at each mesh-point, x = x^_^,y=y 

J 

The equations for computing the exact values of the first horizontal 
derivatives ■ are obtained directly from equation (l-ll) and they are given below: 
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Results and Conclusions 

The exact values of the horizontal derivatives of the field components are 
calculated with the help of equations (l-24) - (1-29). The field components 
are also obtained from total field observations with the method described in 
the previous section. These values of the field components are used to determine 

bicubic spline coefficients in (l-13) are 
calculated. These coefficients yield the values of the horizontal derivatives 
of the computed field components. 

The exact and computed values of the horizontal derivatives for the 
prismatic model in Figure 1 are obtained at equispaced points of a rectangular 
6k X grid system over a horizontal plane. 

Figures 17 and l8 show the exact and computed solutions of respectively. 

C VC 

These two diagrams are almost identical. Some differences between the two 
solutions are noticeable when the magnitude falls off to about one-seventh of 
the maximum response. Fig. 19 shows the K - S and E - ¥ profiles of exact and ' 
computed solutions. The profiles run through the point of maximum response. 

The E - W profile shows the excellent match of exact and computed solutions. 
However, there are appreciable differences between the solutions on the H- S 
profile. At the edges, the computed values on the N - S profiles do not vary 
smoothly. 

Since v V h -0 in the plane of observation, we have 
(1-30) 7) 
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( 1 - 32 ) 

Because of above relations, it is not necessary to calculate all the horizontal 
gradients of the field components. 

Figures 20 and 21 show the exact and computed values respectively of 
Figure 22 shows the .computed solution of *2 Wy /tx- » These three figures are 
almost identical. It is interesting to note the closeness of the computed 
solutions of and This characteristic can also he seen from 

the profiles in Figures 23 and 2h. All the previous diagrams show the discrep- 
ancies between the exact and computed solutions after the magnitudes of the 
derivatives fall off to one-tenth of their maximum amplitudes. The H - S profiles 
contain more discrepancies than the E - W profiles do. For the K - S profile, 
the computed solution of is not as smooth as the computed solution of 

Figures 25 and 26 show the exact and computed solutions of /-q 

respectively. Figure 27 shows the K - S and E - W profiles. All these 
“three .diagrams indicate. that' the values .-are practically identical except in 
regions where the values are less than one-tenth of the maximum amplitude. The 
profiles also indicate that the values of the exact and computed derivatives 
are in better harmony in the case of the E - W profiles than they are in the 
case of the N - S profiles. 

The exact and computed solutions of are shown in Figures 28 and 29 

respectively. The profiles along the N - S direction and E - ¥ direction are, 
shown in Figure 30. Similar quantities are plotted in Figures 31, 32, and 33 
for * Tile results in these diagrams are the same as those in the previous 

figures. A study of all the plots presented here show that of all the derivatives 
the computed has the best match with the exact values. 
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5. Second Horizontal Derivatives of Field Components 

The hicubic spline coefficients, as computed in the previous section for 
determining the horizontal gradients of field components, are also used to cal- 
culate the second horizontal derivatives. Data for the dimensions and magnetiza- 
tion of the prismatic body are given in section 1 of this chapter. 

For the three components of the field, nine second horizontal derivatives 
can be calculated. However, it is known from the derivation in section three 

that • From this relationship we can derive: 
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Therefore, only seven horizontal derivatives have to be calculated. 
However, because of the possibility of inaccuracies in computation, all the 
nine second horizontal derivatives have been calculated. 

The exact values of the second horizontal derivatives are derived directly 

from equation (l-ll) . The equation of these exact values are given below: 
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6. Results and Conclusions 

Figures 3^ and 35 show the exact and the computed values of . 

It -is apparent from Figure 35 that only about one-third of the calculated values 
■which cover the central part of the contoured map, is reliable. However, even 
in this central part, the discrepancies between exact and computed values are 


appreciable. Figure 36 presents the N - S and E - ¥ profiles of the computed 
and exact values. The profiles give a clear picture of how the computed solution 
of ■''3 behaves and why along the boundaries of Figure 35 no contour . was 

drawn. Since the amplitudes of second derivatives fall off sharply, the smallest 
amplitude of reliable data is about one-sixth of the maximum amplitude. Some 
differences can be observed even at the center of the diagram. 

Figures 37 and 38 show the exact and computed values of / 1/ 
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Some differences “between computed and exact fields are noticeable, even in the 

S. 

central part of the diagram. The computed map is not exactly symmetrical. The 
two diagrams show good similarity until the amplitude of the computed field falls 
off to about ' one-fourth of its maximum amplitude. Figure 39 shows the computed 
values of ? . As we pointed out before, ^ is theoretically 

equal to "b . However, the computed values of are not 

exactly equal to those of /'d -x; 'i>y . Furthermore, the computed 

IS closer to the exact values than computed c Hy / "d k . . This characteristic 

is clearly indicated in the N - S and the E - W profiles shown in Figures and 
kl. These two diagrams also indicate the occurrence of the maximum value at 
different points. For the computed solution, the maximum amplitude of computed 

occurs at the point x = -8, y = whereas that of * 

at X = -8, y = 3. 

Figures k2 and present the exact and computed values of 

2 

respectively. The fields of /'0 contain two distinct anomalies, one 

located at the south end and the other at the north end of the prism. The 
computed values are close to the exact values only near these anomalies. These 
anomalies are also present in Figure kk. 

Figures 45 and 46 show the exact and computed values of 
These two diagrams are practically identical except for small values of the 
contours. The N - S and E - W profiles of 0 / ‘^y are presented in • 

Figure 48. As mentioned in section 5, is theoretically equal 

to / 'd X. fl y . Therefore, Figure 45 also represents the exact 

map of 'ey • The computed solution of /%->c '2y is 

shown in Figure 47. Again Figures 45 and 4 t are identical. Figure 49 shows the 
N - S and E - W profiles of 'c Hy / ox . The computed 'i is very 

close to the' exact values. 
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The exact and computed values of y/'Cy are presented in 

Figures 50 and 51. The difference between the two diagrams becomes obvious when 
the amplitude drops to about one-tenth of the maximum. Figure 52 'shows the 
N - S and E - ¥ profiles through the maximum. These profiles also indicate the 
inaccuracies of the computed fields around the boundaries when the amplitude 
becomes relatively small. 

Figures 53 and 5^ show the exact and computed values of ^ ^ 

Some differences between the two diagrams are noticeable when the amplitude 
drops to about one-eighth of the maximum value. The computed field generally 
does not falloff as sharply as the exact field. Figure 55 shows the N - S and 
E - 'W profiles . 

'Figures 56 and 57 present the exact and computed values of • 

The characteristics of these diagrams are very similar to those observed in 
the previous figures. The profiles through the maximum are shown in Figure 58 . 

7 . Concluding Remarks 

In this chapter, the results of investigation on the accuracy of computing 
field components from total field observations have been presented. It is found 
that even for a simple model such as a prismatic body, the computed field com- 
ponents and their derivatives are accurate only in regions surrounding the 
aaaximum amplitudes . As the amplitudes of the various derived fields start to 
decay, the errors in computation increase. 

In the case of actual measurements of the total field, the accuracy in 
computing field components will be much inferior to that obtained in this study. 
The gradients and second horizontal derivatives of the computed components will 
show large errors. These errors will also be responsible for inaccuracies in 
the computation of parameters of causative bodies. It can, therefore, be safely 
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concluded that measurements of field components, if and when possible, is much 
preferable to the calculation of these components, with the help of total field 


observations . 
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Chapter 2: Calculation of the Magnetization Vector 

In the previous chapter, the magnetic field components along three mutually 
orthogonal directions vere obtained from the total field measurements. These 
Components are used to determine the direction-cosines and the magnitudes of the 
magnetization vectors associated vith causative bodies. These unknotm parameters 
are also computed with exact values of the field components due to prismatic 
bodies. A comparison of the values of the parameters, as obtained from exact 
and computed field components, is useful to evaluate the need for measurements 
of field components. 

1. Method for the determination of the magnetization vector 

The method adopted in this study for the determination of the magnetization 
■vector is essentially based on the theoretical treatment of this subject by 
Bhattacharyya ( 1967). 

Let us assume that over the volume occupying the magnetized masses the 
Tnagnetic dipole moment per unit volume has a fixed direction but variable 
magnitude. If the direction is supposed to have an inclination I and declination 
L, and If the direction cosines are denoted by , and V , we have 
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Let us assume that the total magnetic moment of the magnetized mass in volume Vo 
be denoted by 
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C2-2) 
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(2-3) 

(2-4) 

(2-5) 

.( 2 - 6 ) 

(2-T) 

The declination and inclination of magnetization vector can "be easily 
expressed as 

(2-8) y> = ( ^/ v'c<’--+^ ) 

( 2 - 9 ) X r: "/ , 

Thus in the case of three-component data, the magnitude and direction 
of the magnetization vector can he evaluated with the help of equation (2-3) to 
(2-9). 

As a matter of fact, only the moments of H and H in the above equations 
are used here. From equation (2-3) it is easy to see that, we can use the moment 
of H instead of H . Thus it can be concluded that H together with only one of 

y X Z 

the two horizontal components is sufficient for calculating the magnitude and 
direction of the magnetization vector. 
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KoWj let us assume that the total available range of integration in 
the :x direction runs from X to X ^ My^/z and in the y direction, 

from ~ when the origin of the coordinate system 

is placed somewhere near the center of the anomaly, where 
M is the total numher of samples along the x axis, 

N is the total number of samples along the y axis, 

4 is the interval between the samples. 

Therefore, from equation ( 2-3), (2-4) and (2-5) can write 


( 2 - 10 ) 

( 2 - 11 ) 

(2-12) 


Z rr X 
Zn 


.V 


X f’i-y 


•)C 


'X -V . 
X .X 


r>y, 



f-f.. 6Vy 

it 


r M 




-X 

.X 



H, — 


- M. 


-X -Y 


The Integrals ( 2-10) - (2-12) are of the type 

/ / / (■’S r) - 

Using the trapezoidal rule we have 

.Y 


(2-13) 


Y _ ^ 

FOO - / ‘•Y ^ ^ / tYxo 2. Y-,1 z 4 m ' 

^ o' L iir/ -j 


(2-14) 


Then using the same trapezoidal rule for the integrals F(x), we have 

X . _ ^J _ / 
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Assuming that f(x,y) tends smmothly to a very small value along the 
boundary, we can obtain the following approximation from (2-13) and (2-l4): 


X Y /i^l M 
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o-*T t *»TO ^ 


(2-15) 
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vhere f = f Gms, ns) 
jnn 

By assuming unit grid spacing and vith. equalities replacing the approxima- 
tion signs,, the reasonable computational versions of eciuation (2-l0) to (2-12) 


are given by 


M A/ 


(2-16) 
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= ZI 

= I VI ; 

(2-17) 


M ^ 

= 11 
7n= f 

(2-18) 
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With the availability of digital computers, the numerical calculation 


of the above equation is not very time-consuming. 

2, Calculations and Results 

Lourenco Cl9T2) has calculated the magnitude and direction of the 
magnetization vector using exact and computed values of three-component 
magnetic data for different types of theoretical models. 

In the present study also, a rectangular prismatic body is chosen as 
the model and the magnetization of the body is assumed to be arbitrary. The 
parameters selected for the model are given below: 

2a = l6, 2b = 8, h^ = 6, h^ = 10, = i* 

The fo3J-Owing two cases are then studied for examining the influence of the 
geomagnetic field on the calculation of the directions of the magnetization vector: 
(i) I = 60 ° and D = 

(ii) = 30° and = 15®. 


2T 


Por the jaagnetdlzatlon vectorj the inclination is initially fixed at 20° 
and the declination is varied from 0°. to' 80°. at steps of 10°. Then with the 
help of exact and computed field components, the direction and intensity of 
the magnetization vector are determined with the method described in the previous 
section. Studies have also "been conducted for the case when the declination 
is 0° and the inclination changes at steps of 10° from 0° to 80°. Exact and 
computed field components are obtained over a 64 x 6^ grid for the calculation 


of moments . 

The results of calculation are given in tables (2-l) and (2-2). In 
these two tables, D, I and are the assumed values of declination, inclination, 
and intensity of magnetization; and are calculated values for exact 

three-component magnetic data; Dg, Ig and I^g are obtained from computed field 
components for the case, I^= 6o° and = 45°; and D^a 1^ ^p 3 also 

obtained from computed field components for = 30° and = 15°. 

The deviation of declination which is the difference between the 
assumed and computed values of declination, is shown in Figures 59 and 60 as 
functions of inclination and declination respectively of the magnetization vector. 
In Figure 59, the declination has been kept fixed as shown in Table 2-1. Figure 
60 corresponds to Table 2-2. 

Two characteristics of the curves in Figs. 59 and 60 are to be noted. 
First, the computed values of declination deviate with increasing error from 
the correct values as the inclination of the magnetization vector increases. 


Second, errors are significantly larger for lower geomagnetic latitudes. A 
comparison of the two curves for I = 60°, D = 45° and I = 30°,. D = 15° shows 

^ O O O Q 

that the computed declinations vary by about 3-5° over the entire range of 
variation of the direction of the magnetization vector. 
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Figures 6l and 62 present th.e deviation oi* the computed values of inclina- 
tion from the. correct values as functions of inclination and declination respectively. 
In this case also, errors in computation increase as the geomagnetic latitude 
of the observation point is decreased. Errors in values obtained from computed 
field components decrease vith increasing inclination and declination of the 
magnetization vector, become zero and then start to increase vith opposite signs. 

•No explanation is available for this particular characteristic of the errors. 

Tables 2-1 and 2-2 shov that the computed intensities of magnetization 
are practically constant for all possible directions of' the magnetization vector 
and differ by a constant from the correct value. This difference is mainly due 
to the limited area of observations chosen for the calculation- of moments. 

3. Concluding Remarks 

The moments calculated from three-component magnetic data are not 
affected by changes in inclination or declination of the magnetization vector. 

Perhaps moments could have been calculated vitli greater accuracy by choosing 
larger areas with increase in data points. 

With exact values of field components the error in calculated declinations 
is always less than two degrees and the computed inclinations are within 1° of 
the exact values. The computed field components, as expected, give rise to 
greater errors and these errors increase with reduction in geomagnetic latitude. 

The calculations presented in this report show rather conclusively the 
superiority of exact components over computed components in determining the 
magnetization vector with a reasonable accuracy. 


Table 2-1 


D = 0°, 1=1 
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I 


^1 ■ 

^pl 

“2 

I 2 

^P2 

^3 


■ I -3 
p 3 

0 

' 0.002 

0.795 

0.665 

0.80 

5-211 

0,645 

4.153 

7.391 

0.561 

10 

0.083 

10.786 

0.666 

1.194 

14,058 

0.658 

4.618 . 

18.178 

0.584 

20 

0 . 1 T 3 

20.760 

0.666 

1.606 

22.629 

0.665 

5.100 

28.112 

0.609 

30 

0 . 2 T 6 

30.720 

0.667 

2.067 

31.079 

0 . 668 . 

5.633 

37.297 

0.632 

ho 

0 . 1+03 

40.672 

0.667 

2.619 

39.552 . 

0.663 

5.265 ■ 

45.888 

0,650 

50 

0.576 

50.620 

0.667 

3.344 

43.194 

0.654 

7.080 

54.063, 

0.663 

60 

’ 0.841 

60.571 

0.667 

4 . 4 ii 

57.143 

0.640 

8.256 . 

62.001 

0.667 

TO 

1.347 

70.530 

• 0.666 

6.284 

66.519 

0.623 

10.251 

. 69.875 

0.658 

80 

2.856 

80.498 

0.666 

, 10.853 

76.383 

0.605 

14.769 

77.830 

0.667 
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D 

^1 

^1 

^pi 

^2 

"a 


"3 

S 

^P 3 

0 

0.173 

20.706 

0.666 

1.606 

22.629 

0.665 

5.100 

28.112 

0.609 

10 

9.962 

20.793 

0.-666 

11.494 

22.184 

0.668 

16.149 

27.341 

0.619 

20 

20.090 ' 

20.790 

0.667 

21.265 

21.627 

0.670 

26.639 

26.306 

0.632 

30 

30.197 

20.754 

0.667 

30.933 

20.983 

0.670 

36.529 

25.119 

0.645 

1+0 

40.374 

20.688 

0.668 

40.518 

20.274 

0.670 

45.857 

23.878 

0.656 

50 

50.317 

20.596 

0.669 

50.051 

19 - 524 

0.668 

54.714 

22.658 

0.665 

60 

60 ..324 

20.486 

0.670 

59.566 

18.752 

.0.665 

63.223 

21 . 511 

0.670 

70 

70.299 

20.362 

0.670 

69.100 

17.978 

0.659 

71.520 

20.470 

0.670 

80 

80.247 

20.232 

0.670 

78.689 

17.222 

0.654 

79 . 7 i ^8 

19.555 

0.665 

90 

89.982 

20.100 

0.670 

88.365 

16.493 

0.647 

88.055 

18.777 

0.654 
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